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Abstract 

In the framework of Landau theory of phase transitions one is inter- 
ested to describe all the possible low symmetry "superconducting" phases 
allowed for a given superconductor crystal and to determine the conditions 
under which this crystal undergoes a phase transition. These problems 
are best described and analyzed in the orbit space of the high symmetry 
group of the "normal, non-superconducting" phase of the crystal. In this 
article it is worked out a simple example concerning superconductivity, 
that shows the P-matrix method to determine the equations and inequal- 
ities defining the orbit space and its stratification. This approach is of 
general validity and can be used in all physical problems that make use 
of invariant functions, as long as the symmetry group is compact. 1 

1 Introduction 

Landau theory of phase transitions 15] Hj has been sometimes used to 
study superconductor crystals. In that framework one deals with a 

given potential function V(x), defined in a finite dimensional vector space 
X, that is invariant under transformations of a compact group G acting 
linearly on X: V(gx) = V(x), V<? £ G, x £ X. The vectors x £ X are 
called order parameters. All finite dimensional representations of a com- 
pact group G are unitary and completely reducible, and, given a complex 
irreducible representation <p of G, either ip itself, or the direct sum ip + ip, 
where (p is the complex conjugate representation of tp, is equivalent to 
a real orthogonal irreducible (on the reals) representation of G. |S] Q I n 
Landau theory, in a first approximation, sufficient to determine all possi- 
ble symmetry phases allowed for the system, one considers X hosting a 
representation of G equivalent to a real orthogonal irreducible representa- 
tion. 2 This implies that no vector x is left invariant by all elements 
of G, so the action is effective. In all generality, one may then assume that 
X = R n and that the group G is a real orthogonal subgroup of 0(n, R). 
In the following the irreducibility condition will never be used, but we still 
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require effective orthogonal actions in R™. 

In this Introduction some of the mathematical tools developed in are 
reviewed; they are used in Section 2 to determine the symmetry phases 
allowed for the physical system considered. The method exposed modi- 
fies slightly that one proposed by Gufan. |1U| The physical results here 
obtained are not new, but I would instead point the attention of the 
reader to the mathematical method used to derive the orbit space strati- 
fication and the list of phase transitions allowed for a given potential. 

If xo is the minimum point of V(x), then the system is in equilibrium 
in a state compatible with the order parameters xq. If xo 7^ then the 
symmetry of the system is lower than G, it is reduced to G xo — {g G G | 
gxo — xo}, the isotropy subgroup of xo, realizing a spontaneous symmetry 
breaking. 

Moreover, if xo 7^ 0, the orbit Q(xo) = {gxo, g G G} contains points x 7^ 
xo where V(x) — V(xo), so the low symmetry equilibrium states are al- 
ways degenerate. The following relation holds: ^JE] G gx = gG x g~ , Vg G 
G, x £ R n , then all points x' G Q(x) have isotropy subgroups in the same 
conjugacy class [G x ] and viceversa, all groups G' G [G x ] are isotropy sub- 
groups of some points x G Q(x). The conjugacy class [G x ] is called the 
(isotropy) type of Q(x). Two conjugated isotropy subgroups are physi- 
cally equivalent because they can be thought to be the same group seen 
from two rotated systems of reference, of course we are here considering 
G orthogonal. Given an isotropy subgroup H C G, the set of points 
Erm = {x G R n I G x G [H]} is called the stratum of type [H]. The strata 
are disjoint and in a one to one correspondence with the symmetry phases 
that may in principle be accessible to the system. Orbits and strata can 
be partially ordered according to their types: the type [H] is said to be 
smaller than the type [K], [H] < [K\, if H' C K' for some H' G [H] and 
K' G [K]. If G is compact and is acting orthogonally in R n , the number 
of strata is finite, | I3| p. 444, and there is a unique stratum of the smallest 
type, |12| called the principal stratum. 

The values assumed by the potential V generally depend, in addition to 
x, on other variables, like temperature, pressure, atomic concentrations, 
etc., and the absolute minimum of V falls in a point xo that depends 
on these variables. Sometimes, when these variables change their values, 
xo may change stratum, and in this case the system undergoes a phase 
transition. 

The orbit space of G is the quotient space R n /G, whose points rep- 
resent whole orbits of R". A stratum of the orbit space is formed by all 
the points that represent orbits in a same stratum of R". The invariant 
functions, e.g. V(x), axe naturally defined in the orbit space because they 
are constant on the orbits, so, to eliminate the degeneracy associated to 
the symmetry, one is lead to study the invariant functions in the orbit 
space. 

To determine concretely the equations and inequalities that define the 
orbit space and its strata one may take advantage of the P-matrix and 
its properties. HUE] By Hilbert's theorem, P3[T?>] the algebra R[R"] G 
of real G-invariant polynomials defined in R n has a finite number q of 
generators: pi, . . . ,p q G R[R"] G , forming an integrity basis for R[R n ] G . 
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This means that 

f(x) = f( Pl (x),..., Pq (x)) V/£R[Rf, x € R™, 

with / a real polynomial in q variables. Actually, the preceding formula 
holds true for any G-invariant C°° function f(x), then / is a C°° function 
of q variables. |17l I18 | If no subset of {pi,...,p q } is still an integrity 
basis for R[R"] G , then the set {pi, . . . ,p q } forms a minimal integrity basis 
(MIB) for R[R"] G . There is some arbitrariness in the choice of the MIB, 
however all its elements can be chosen to be homogeneous and one may 
order the p a , a = 1, . . . ,q, according to their degrees d a = Aegp a (x), for 
example d a > d a +i- If G acts effectively in R™, then d q — 2, and, for the 
orthogonality of G, one may take p q (x) — \x\ 2 = ^™ =1 x i- 

The vector map p : R" — » R 9 : x —* p(x) = (pi(x),p2(x), . . . ,p q (x)), 
called the orbit map, is constant in each orbit Q C R" because the p a (x) 
are G-invariant functions. The point p = p(x) G R 9 is the image in R 9 of 
the orbit Q(x) of R". No other orbit Q' 7^ Q of R n is mapped to the same 
point p £ R 3 because the MIB separates the orbits. 18 8 The image of 
R" through the orbit map is the set 

S = p(R") C R 9 

that can be identified with the orbit space of the G-action. |18| 
The strata a of S are the images of the strata E of R n through the orbit 
map: a = p(E). The principal stratum a p of <S is connected, open and 
dense in S. |19| Given two strata E and E' of R n , if E' is of greater type 
than E, then a' — p(E') lie in the boundary of a = p(E). |20| The smallest 
stratum of 5 is <t\q\ , image of the origin of R n and is located at the origin 
of R 9 . Given any plane n r of R 9 with equation p q = r 2 > 0, 5 n Il r , is 
a non-empty compact connected section of S that contains all strata of S 
except <7r<3] . |5] P| This section is sufficient to describe the whole shape of 
S, because moving Il r toward the infinite or toward the origin of R 9 , by 
changing r 2 , 5 n H r expands or contracts, but maintains its topological 
shape and stratification. 

Let's define the qxq symmetric and positive semi-definite matrix P(x), 
with elements 

D . . , . , . dpa(x) dp b (x) 

Pab(x) = X7pa{x) ■ Vpb(x) = —q^ q^— a,b=l,...,q. 

i=l * 1 

The matrix elements P a b(x) are real homogeneous polynomial functions of 
x with deg P a b(x) = d a + d b — 2, and P qa {x) = P aq (x) = 2d a p a (x). A less 
immediate property is the G-invariance; it follows from the orthogonality 
of G and from the covariance of the gradients of G-invariant functions: 

Pab{gx) = Vpa(gx)-Vp b (gx) = gX7pa(x)-gX7p b (x) = \7pa(x)-Vp b (x) = P a b{x). 

Then, P a b(x) can be expressed in terms of the MIB: 

Pab(x) = P a b(pi(x), . . . ,p q (x)), Vi 6 R°, a, b = 1, . . . , q. 

The matrix P(p), defined in R 9 , is called the P-matrix. When p — p(x), 
then P(p) coincides with P{x). 
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Let's call a polynomial f(p) w -homogeneous of weight d if the poly- 
nomial f(x) — f(p(x)) is homogeneous of degree d, and define the sur- 
face of the relations: Z — {p £ R 9 j fi(p) = 0, i = 1, ...,&}, where 
fi(pi(x), . . . ,p q (x)) = 0, i = l,...,fc, are fc algebraically independent 
homogeneous polynomial relations, called syzygies, existing among the 
p a (x). If there are no syzygies, then k = and Z = R 9 . 

The main properties of a P-matrix are the following: |H| 12 "I | 

• P(p) is a real, symmetric q x q matrix. The matrix elements P a b(p) 
are w-homogeneous polynomials of the p a of weight d a + db — 2, and 

Pqa(p) = Paq(p) = 2d a p a , Vtt = 1, . . . , £/. 

• P(p) is positive semi-definite in S. Precisely, S = {p G Z | P(p) > 

o}; 

• Given a stratum a C S, and p £ <r, dimcr = rankP(p). 

The P-matrix completely determines the orbit space S and its stratifi- 
cation. Defining Sd as the union of all d-dimensional strata of <S, one 
has: 

S d = {p G Z | P(p) > 0, rankP(p) = d}, 

so, <S = S q -k, where q — k = dim Z and k is the number of syzygies. 
As all these positivity and rank conditions are expressed through a finite 
number of polynomial equations and inequalities, S is a semi-algebraic set 
with a finite number of strata (confirming Mostow result |13p. 

The polynomials defining the strata a of S of dimension q — 1 (that 
exist only if there is at most one syzygy) must satisfy a set of differential 
relations that characterize them. 1221 Let a(p) be an irreducible poly- 
nomial defining such a stratum of S. Then the following master relations 
hold true: 

y^Pabip)^^- = \ a (p)a(p) a = l,...,q, 
opb 

b=i 

where the A a (p) are to-homogeneous polynomials of weight d a — 2. The 
polynomials a(p) that satisfy the above master relations are called active. 
They are always factors of detP(p). |3] The (q — l)-dimensional stratum 
a is located in the interior of a — {p G R 9 | a(p) = 0, P(p) > 0}. The 
border a\a of a is the union of lower dimensional strata whose defining 
equations and inequalities can be found simplifying the system {ai(p) = 
£i2 (p) = 0, P(p) > 0}, where oi(p) and a2(p) are two irreducible active 
polynomials. 

As the P-matrices completely define the orbit spaces and their strati- 
fication, the classification of P-matrices implies the classification of orbit 
spaces. All P-matrices of dimension q < 4 corresponding to coregular 
groups (i.e. with no syzygies) have been determined and classified. |51l23| 
This classification is obtained independently from the specific group struc- 
tures by making use only of some general properties of orbit spaces. It is 
found that many different groups, no matter if they are finite groups 
or compact Lie groups, have the same P-matrix and the same orbit 
space. |2U 
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2 A concrete example analyzed in the or- 
bit space 

The example here proposed was first studied by Gufan and concerns 
the study of the possible symmetry phases of electron pairs with non-zero 
total orbital momentum in hexagonal superconducting crystals, in the 
approximation of strong crystalline field and strong spin-orbit interaction. 
In |1 1 1 the order parameters are chosen to be two complex numbers, r/+ 
and r\- , and the symmetry group is Y — G ® R Cg> Ui(a), where G is the 
point symmetry group of the crystal, 7? is the time reversal group and 
Ui(a) is the rotation group of the complex plane. Here we consider only 
the case G — C?, v . The group Y acts in the four dimensional complex 
space of vectors of the following form: 




where * denotes complex conjugation, and the group generators have the 
following matrix expressions: 








where C\ z and Ui(a) are diagonal matrices, of which only the diagonal 
elements are shown. 

The MIB for this group has 3 elements that can be written in the following 
real form: 

J:i = \v+\ 2 + \V-\ 2 , J2 = \v+\ 2 \V-\ 2 , Ji =V+V- 3 + V+ 3 V-- 

Of course one may consider the action of Y in C 2 , but this action would 
not be linear, the invariants would not be polynomial functions of rj + and 
r\- , and the theory exposed in the Introduction would not be valid. 
When a unitary representation of a compact group is equivalent to a real 
orthogonal representation, one may always choose the MIB's for the two 
cases (unitary and orthogonal) to be formed by the same real polynomial 
functions. |24| This means that the orbit map defines in the two cases the 
same subset S C R 9 , that is the same orbit space. 

To adapt this scheme to that one reviewed in the Introduction, one has 
to consider an orthogonal action of Y in the real space R 4 . Introducing 4 
real variables £1,2 and 2/1,2, one may put rj+ — x\ +iyi and r)- = £2 + *J/2, 
and consider a vector x € R 4 of the following form: 
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This defines a one-to-one correspondence between the complex space of the 
77 and ft 4 . The group Y acts in R 4 with the following matrix generators: 



(0 1 \ 

1 I 

10 

10 0/ 

R 

1 \ 
0-11 

0' 
-10 0/ 

The MIB can be easily written down in terms of the real variables: 

p 3 {x) = x 1 + x 2 + y 1 + y 2 , 
p 2 (x) = i(xf +yt)(xl +yl), 

Pi(x) = 16(x!X 2 +yiy2){x(x? 2 ~3xly 2 +8x!X 2 yiy2 ~ Sxlyl + yly 2 ), 

where the numeric factors have been introduced for later convenience. 
To find out the P-matrix one has to calculate all gradients of the p a {x), 
make their scalar products, and express these products in terms of the 
p a (x). One finds, at the end, the following P-matrix: 

/ 144p|p 3 24pip 3 12pi \ 
P(p) = 24pip 3 16p 2 p3 8p 2 , 
\ 12 Pl 8p 2 4p 3 / 

whose determinant is 

det P(p) = 2304 p 3 {p\ - ipl)(p 2 -p 2 3 ). 

The p a (x) are algebraically independent, so the principal stratum is the 
subset of ft 3 where P(p) > 0, and all other strata are contained in the sur- 
face of equation A(p) = 0, where A(p) is the product of all the irreducible 
active factors of det P(p). Using the master relations, one finds: 

A{p) = {p\-Apl)(p 2 -pl). 

The section of the surface A(p) = in a plane P3 =constant is plotted 
in Figure Q Its interior hosts the section of the principal strata a p and 
its border contains the sections of all other strata, except the origin. The 
complete list of the strata of S is reported in Table Q where the defining 
equations and inequalities of the bordering strata are found by examining 
rankP(p) in the surface A(p) — 0. 

Table |21 lists for each stratum a, some conveniently chosen points x and 
the generators of their isotropy subgroups G-. One may note that the 
points in the border of a stratum a have isotropy subgroups containing 
those of points of a, modulo a conjugation. Just to see this fact, in two 
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where c = cos a and s = sin a. 
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Fi gure 1* Section of the orbit space. 



Table 1: Strata of the orbit space 



Stratum 


Denning relations 


CO 


Pl = P2 = P3 = 


<J\ 


Pi = P2 — 0, p 3 > 


01+ 


Pi = 2pl, P2 = p|, P3 > 


Cl- 


Pl = -2p|, p 2 =p§, p 3 > 


02 


|pi| < 2pl, p 2 =Ps, P3 > 


02+ 


pi = 2p|, <p 2 <p|, p 3 > 


CT2- 


pi = -2pj, 0<p 2 <pi, p 3 >0 


a v 


pf < 4p3, < p 2 < pl p 3 > 



cases are reported two different points in a same orbit, related by the 
group element Ui(^). 

Once the orbit space is described, one may list all second order phase 
transitions that may in principle be allowed for the system, a second order 
phase transition may in fact take place only between bordering strata 
because it is a consequence of a continuous variation of the variables that 
appear in the potential. ^ |5J |2| In the example here considered, the 
possible second order phase transitions are the following ones: 

o"o,cTp < > all other strata 02+ < > d, <7\+ 

(72 < > <Tl+, <Tl- <T2- < > CTl, (Tl_ 

The conditions under which the phase transitions are realized depend on 
the form assumed for the potential when this is written in terms of the 
MIB and not all the transitions listed above do actually take place when 
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Table 2: Isotropy types of the strata 



Stratum 


Typical point x 


Generators of the isotropy subgroup of x 




(0,0,0,0) 


U^a); C^; a v ; R 




(0,1,0,0) 


L7i(-§7r)Cg*; U^)a v R 


01+ 


(0,1,0,1) 


ct v ;U 1 {-k)R 


0\- 


(-1,0,1,0) 


Ui(n)a v ; U X {%)R 




(0,1,0,-1) 


Ui{ir)a v ; R 


0"2 


(-1,1,1,1) 


Ui{tt)R 




(1,1,1,-1) 


R 


0"2+ 


(0,1,0,2) 


Ui{-k)(t v R 


02- 


(0,1,0,-2) 


Ui(ir)a v R 


(Jp 


(-1,2,1,1) 


E 



a given form of the potential is specified. It is not possible to examine 
here in detail the phase transitions scenario associated to a given potential 
V(p) and I shall here only remind the basic points of this analysis. |24llll| 
One first has to write down V(p) in terms of the MIB. Generally, one uses 
a polynomial approximation for V(p): V(p) = clips + aip\ + . . . + &1P2 + 
. . . + api + . . ., truncated to a certain degree. One determines then the 
conditions on the parameters a,i,bi,... in order that V(p) be bounded 
below in each of the strata of S, determines these constrained minima, 
and compares them to single out the absolute minimum. At the end one 
knows, for any set of values for the parameters, which is the stable phase, 
that is the stratum hosting the absolute minimum of V(p). It is evident 
that polynomial potentials of small degrees cannot develop minima in 
certain strata of S for any set of values of the parameters. 

The order parameters are linked to some physical observables. This 
is the case, for instance, of the number of pairs n — 77+ 1 2 + |?7-| 2 = f>3, 
of the number z = \r)+\ 2 — \rj~\ 2 = ±v^y>§ — pi that changes sign under 
time reversal and determines the magnetic properties of the crystal, [5] 
and of the components of the deformation tensor of the unit cell of the 
crystal. From Table^one sees that in the strata ao, o\- , <ti+ , 02, and only 
in them, z vanishes identically, so these strata cannot exhibit magnetic 
properties. The change of the sign of p\ does not modify the P-matrix 
and the orbit space S, but changes the strata <7i_ <-> 01+ and 02- 02+ • 
As this sign change has no physical meaning, the strata <ti_ and 02- are 
not physically distinguishable from <7i+ and (72+, nevertheless they are 
distinct strata lying in different positions of the orbit space, and phase 
transitions concerning each one of these strata are possible. 

The P-matrix and the orbit space here described, also appears in the 
case of superconducting cubic crystals, if the order parameters are as- 
sumed to transform according to the two-dimensional complex represen- 
tation of the cubic symmetry group. 125) In that case however the isotropy 
groups are no longer those reported in Table [5] As already noted, many 
different symmetry groups may give rise to the same P-matrix and the 
same orbit space. When this happens, the number of different symmetry 



8 



phases and the list of the phase transitions allowed for these systems are 
the same (but the corresponding isotropy groups are generally different, 
depending on the specific high symmetry group considered). Moreover, 
for a given form of the potential in terms of the M1B, the conditions 
to realize the listed phase transitions impose the same conditions on the 
parameters of the potential. 
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